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ABSTRACT 

In this paper we present a fully relativistic approach to modelling both the continuum 
emission and the reflected fluorescent iron line from a primary X-ray source near a 
Kerr black hole. The X-ray source is located above an accretion disc orbiting around 
the black hole. The source is assumed to be a static point source located on an ar- 
bitrary position above the disc, on or off the axis of rotation. We carry out Monte 
Carlo simulations in order to estimate the iron line spectrum as well as its equivalent 
width. Because of the gravitational lensing effect, an enhancement of the iron line 
is expected when the primary source is located close to the central black hole. We 
find that for a source located on the axis of rotation the enhancement is relatively 
modest. An observer at inclination 30 degrees would measure an equivalent width of 
^300 eV in the extreme case of a maximally rotating black hole and a source located 
at height 1.5 gravitational radius from the centre. This corresponds to an equivalent 
width enhancement factor of about 2 compared to the classical value where no lensing 
effect comes into play. However, when allowing the source to be located off the axis of 
rotation, much stronger enhancement can be obtained. In the extreme case of a max- 
imally rotating black hole and a source located just above the approaching side of the 
disc, an observer at inclination 30 degrees could measure an equivalent width as high 
as ~1.5 keV (i.e. ~10 times the classical value). We also find that observers located at 
high inclination angles observe a stronger line than observers at low inclination angles. 
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1 INTRODUCTION 

It is now commonly agreed that Active Galactic Nuclei 
(AGN) are composed of a central supermassive black hole 
closely surrounded by a thin equatorial accretion disc, the 
emission process being controlled by viscous transport of 
angular momentum (see e.g. Rees 1984; Blandford and Rees 
1992). The disc material is relatively cold (~ 10 B K) com- 
pared to the hard X-ray corona (~ 10 s K) by which it is 
surrounded. The corona is responsible for the observed X- 
ray continuum emission which follows a power-law of index 
r = —1.7. X-ray observation of Seyfert-1 galaxies revealed 
that there exist spectral features and other deviations from 
a simple power-law (e.g. Pounds et al. 1989; Nandra et al. 
1989; Matsuoka et al. 1990). In particular a strong iron line 
is observed at 6.4 keV together with an excess in the contin- 
uum emission at energies higher than ~4 keV. This is well 
explained by Compton reflection re-processing of the hard 
X-ray power-law continuum onto the accretion disc (see e.g 
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George & Fabian 1991; Matt, Perola & Piro 1991 ). The ob- 
served iron line is the result of fluorescent processes and the 
excess in the continuum is accounted for by the presence of 
a continuum component in the reflected emission. The de- 
tailed spectrum of the iron line is of great importance since 
it allows one to probe regions of the accretion disc as close 
as a few Schwarzschild radii from the central black hole (see 
e.g. Tanaka et al. 1995). In the particular case of MCG-6-30- 
15, the lineshape yields evidence that the emission occurs as 
close as half a Schwarzschild radius, giving strong support 
for the presence of a rotating Kerr black hole (Tanaka et al. 
1995, Dabrowski et al. 1997). 

The fluorescent iron line profile is well predicted by 
models which assume a power-law e(r) oc r~ a for the fluores- 
cent emissivity on the accretion disc, where a ~ 2 (see e.g. 
Tanaka et al. 1995; Bromley, Chen & Miller 1997; Dabrowski 
et al. 1997; Reynolds & Begelman 1997; Cadez et al. 1998). 
In such cases the hard X-ray corona is assumed to be an 
extended source located above and below the disc so that 
approximately half of the emitted power illuminates the disc 
and is responsible for the reflected continuum component as 
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well as observed fluorescent lines (George & Fabian 1991). 
The remaining half escapes to infinity and constitutes the 
major part of the observed continuum spectra. The pre- 
dicted iron line equivalent width (EW) is in this case es- 
timated to be ~100 - 200 eV (e.g. George & Fabian 1991, 
Matt et al. 1991). However, as discussed in Nandra et al. 
(1997), for a sample of 18 Seyfert-1 galaxies, observations 
reveal equivalent widths of 300-600 eV, or even as high as 
1 keV in the case of MCG-6-30-15. Part of the EW enhance- 
ment may be explained by an overabundance of iron in the 
accretion disc material (e.g. Lee et al. 1998, 1999). How- 
ever, other means of enhancement need to be found in order 
to account for large values of the iron line EW (Reynolds 
& Fabian 1997). Following the work of Martocchia & Matt 
(1996) and very recently Martocchia, Karas & Matt (1999), 
in this paper we propose to investigate the possible EW en- 
hancement due to the gravitational lensing of the primary 
light rays which illuminate the disc and therefore drive the 
fluorescent process. The primary hard X-ray source is as- 
sumed to be a static point source located above the accre- 
tion disc. Situations where the source is located on or off the 
axis of rotation are both considered here. 

Our predictions are based upon Monte Carlo simula- 
tions which account for the primary X-ray emission as well 
the reflected iron line in a consistent manner. It is assumed 
here that the primary continuum emission dominates over 
the reflected continuum emission so that the reflected con- 
tinuum is neglected while calculating the line EW (see Sec- 
tion ^]). Details of our fully relativistic theoretical approach 
are given in Section |^. Some results are given in the form 
of iron line profiles and equivalent widths in Section ^ for a 
primary source located on the axis of rotation. In Section 
this assumption is relaxed and we investigate line profiles 
and equivalent widths for arbitrary positions of the source 
above the disc. Finally our results are discussed in Section ^| 



2 REFLECTION MODEL IN THE KERR 
METRIC 

The black holes considered here belong to the Kerr family 
of solutions to Einstein's equation and the metric employed 
is in the Boyer-Lindquist form, which is defined by 

i ,,i j / dr 2 ■_, \ , , i \ i /, i ■_' 



ds z = dr-p" ( -j-+d9 }- (r' + L*) sin 



2Mr 



L sin 6d<j) — dtj 



where 

p 2 = r 2 + L 2 cos 2 8, 
A = r 2 -2Mr + L 2 . 



(1) 

(2) 
(3) 



Here t, r, 8 and 4> are the space-time coordinates, M is the 
geometric mass, the mass of the black hole being Mc 2 /G 
and the quantity LM is the angular momentum of the hole 
as measured at infinity (e.g. D'Inverno 1992). The gravity of 
the disc itself is assumed to be negligible and therefore the 
whole spacetime is described by the Boyer-Lindquist metric. 
We employ natural units G = c = % = 1, unless stated 
otherwise. 



2.1 Photon Path 

Throughout this paper we use the gauge-theoretic formal- 
ism of Lasenby, Doran & Gull (1998) which allows fully rel- 
ativistic and covariant calculations to be made in the con- 
text of an easy-to-handle flat space. It is important to note 
that predictions for measurable quantities agree exactly with 
those of General Relativity for this case. Equations relevant 
to an equivalent approach to the problem, but following a 
more standard general relativistic formalism, are presented 
in Appendix 

The gauge-theoretic approach to gravity of Lasenby et 
al. (1998) employs the language of geometric algebra, which 
seems to most clearly expose the physics involved. In par- 
ticular, the geometric algebra of spacetime (Hestenes 1966), 
known as the Space Time Algebra (STA), is generated by a 
set of four orthonormal vectors {7,1}, p. = 0...3, satisfying 



7m -7" = V»v = diag(+ ) 



(4) 



In order to investigate problems in relation with the 
Kerr geometry, we also introduce a suitable set of coordi- 
nates t, r, 8 and 4> associated with the polar coordinate 
system and defined in terms of the fixed {7 M } frame. The 
corresponding coordinate frame e M = d^x is 

et = 70 (5) 

e r = sin 8 cos 1^71 + sin 8 sin ^72 + cos #73 (6) 

eg = r (cos 8 cos 07i + cos 8 sin ^72 — sin 673) (7) 

= rsin^ (— sin07i + COS072) . (8) 

Note that eg and e$ are not unit vectors and differ by factors 
r and rsin# from the usual spherical-polar basis vectors, 
respectively. This is a direct consequence of their definition 
in terms of the {7 M } frame (see Lasenby et al. 1998). 

The translation and rotation gauge fields correspond- 
ing to the Boyer-Lindquist metric (|l]) are given in Doran, 
Lasenby & Gull (1996). 

We start by defining an observer with covariant 4- 
velocity equal to 70. Hereafter this observer will be called the 
70-observer. One can show that, at infinity, the 70-observer 
is in a flat Minkowski space-time and at rest with respect to 
the central black hole. Elsewhere its function is to provide 
a useful frame in which to express quantities. In order to 
parameterise the energy and direction of the photon path, 
let us define the photon 4-momentum p in the {7,1} frame 
by 



p = 



$70 + $ sin 8 P cos 4>vli 
<E>cos0p73. 



+ $ sin 8 P sin ^72 + 



(9) 



The photon energy as measured by the 70-observer is given 
by p ■ 70 which, from equation (|^), is equal to $. Similarly, 
as illustrated in Fig. |l| the angles 8 P and (f> p correspond to 
the usual spherical-polar angles in the local spatial frame of 
the 70-observer and define the direction at which photons 
are received (or emitted) by the 70-observer. We note that 
p is guaranteed to be null since p ■ p = 0, which is what is 
required for a massless particle. The photon trajectory itself 
is parameterised by the affine parameter A. The spacetime 
position of the photon (t, r, 8, (f>) is defined by the position 
vector 



te t + re r 



(10) 
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/ Figure 2. The 4- velocity of the orbiting particle, as measured by 

the 70-observer at position (r, 9 = tt/2, <fi). lot = n/2 + <f> and 
Figure 1. Photon 4-momentum p as measured by the 70- F 70 = tanh C/ are respectively the measured angle and magnitude 

observer. 9 P and <p p are the observed angles of the photon tra- Q f ^ ne ve ] oc ity 

jectory. 



where t, r, 9 and cj> are scalar functions of A. The photon 
geodesic equations are given in Appendix ^. The equations 
we obtain are all first order ordinary differential equations 
for t, r, 9, <f>, <E>, 8 P and <f> p . As mentioned in the appendix, 
the attractive aspect of this approach is that the geodesic 
equations have been obtained without making use of con- 
served quantities along the photon path such as the energy, 
angular momentum or Carter constant. Instead these quan- 
tities can be used a posteriori as an effective check of both 
analytical and numerical results (see Appendix |§|) . 

2.2 Accretion Disc 



We now employ a similar method as in Section 2.1 to study 



the properties of an accretion disc rotating around a Kerr 
black hole. In particular we are interested in the velocity 
and energy of the accreted material. The assumptions are 
the following: 

(i) The accretion disc is a thin disc and lies in the equa- 
torial plane (9 — tt/2), perpendicular to the axis of rotation 
of the central Kerr black hole. 

(ii) The orbits are stable circular orbits. In the same 
way as we defined the photon momentum Q), we start here 
by defining the 4- velocity vector Vd of a particle in circular 
orbit of radius r around the black hole. Since the motion 
occurs in the equatorial plane, we can parameterise Vd as a 
function of 70, 71 and 72 only: 



v d 



cosh (770 + (sinh (7 cosojt)7i + (sinh (7 sintjr)72, (11) 



where u) represents the angular velocity and r is the proper 
time along the world-line of the particle in circular orbit. U 
is a function of r, L and M, but not a function of r since 
the orbit is circular. The parameterisation in cosh and sinh 
has been chosen to ensure that 



Vd - rid = cosh 2 U — sinh 2 (7 = 1, 



(12) 



as required for a massive particle. In the {e M } frame, the 
velocity takes the simple form: 



Vd = cosh Uet + 



sinh U 



(13) 



We note that the 70-observer is not comoving with the 
particle in orbit. Let us however assume that, at a given 
time t, the 70-observer is at the position (r, 9 = n/2, <f>), i.e. 



instantaneously at the same position as the orbiting particle 
(see Fig. fy. According to the 70-observer, the energy per 
unit mass of the particle is given by E 10 = Vd • 70 = cosh U. 
The particle velocity has a magnitude Vy = tanh(7 and 
makes an angle lut — n/2 + <f) from the 71 axis in the local 
spatial frame (71, 72, 73) of the 70-observer. When solving 
the dynamical equations, we find 



tanht/ = 



-L± 



A2 



(14) 



and therefore 



cosh U 



At 



yr 2 



sinh U — 



3Mr ± 2L\ 
-L±yfMr 



=• (15) 

\/r 2 - 3A/r ± 2L\/Jfr 

In the case where L > 0, the + sign in (|l4j) and ([[5]) is for 
orbits co-rotating with the spinning hole, while the — sign 
is for counter-rotating orbits. If L < the significance of the 
signs is reversed. 

Finally, following the same approach as in Appendix [b], 
the conserved energy per unit mass E r and angular momen- 
tum J r of a circular orbit can be evaluated by using the 
Killing vectors Kt (equation |b^) and iQ, (equation |bJ). We 
find 



and 



KfV d = E 1 



r 2 - 2Mr ± L\ 



Jr = K<f, ■ V d = E n 



LA 2 +L 2 V 10 



+ E- ln V- ln r. 



(16) 



(17) 



The conserved quantity E r in ( |l6| ) is related to the local 
energy E 10 by a redshift factor, which tends to unity at 
infinity. Therefore E r is usually called the energy at infinity. 
Regarding J r in (p"7|), the first term is a purely gravitational 
effect due to the dragging effect of the spinning hole, while 
the second term is more directly related to the Newtonian 
expression, mVr, of the angular momentum. We note that 
the condition for stability is given by dE r /dr > 0, which 
gives the following condition for the minimum stable orbit 



tanht/(r ms ) = -. 



(18) 
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Figure 3. Geometrical set up of our Monte Carlo simulations. 
The spinning black hole is located at the centre of coordinates 
and the accretion disc lies on the equatorial plane from the inner 
radius r; n to the outer radius r on t ■ The arrow indicates the sense 
of rotation. The primary source is situated above the disc at an 
arbitrary position (p s , (ps, z s ) and illuminates both the disc and 
the collecting area of a distant observer located at (r , O , <j> ). 
The surface of the collecting area is r% sin O A</> A0 O . 



It is perhaps somewhat surprising that the minimum stable 
orbit is for a local velocity of one half the speed of light. In 
this work, the accretion disc extends from r = r- m — r ms up 
to the arbitrary large outer radius r — r out = 5000 r g , where 
r g = GM/c? is the gravitational radius. 

2.3 Primary Source and Distant Observer 

The primary X-ray source is assumed to be a static point 
source and, as shown in Fig. ^, is located above the disc 
at an arbitrary position (p 3 , <f> S) z s ), where p, <j> and z are 
the usual cylindrical polar coordinates. The 4-velocity v 3 of 
the source has to be proportional to the Killing vector K t 
(equation B3) so that the source can remain static (see e.g. 
Misner et al. 1973). Requiring that v s ■ v 3 = 1, we find 

1 



--K t 



(19) 



The source is emitting isotropically in its proper frame with 
a rate of emission defined by 

N l/s oc z/f photon s" 1 Hz" 1 , (20) 

where V is the photon power-law index. Throughout this 
work we assume T = —1.7 (Mushotzky 1982; Turner & 
Pounds 1989). 

Predictions regarding the iron line flux and equivalent 
width are those measured by a distant static observer. Since 
Kt ~ 70 for large r, the 4-velocity of the distant observer is 

^=70- (21) 

As illustrated in Fig. ^, the position of the distant observer 
is given by r , 9 , <f> with collecting area on the r — r 
sphere equal to 



AA = r a sin 6 A4> A8 . 



(22) 



2.4 Monte Carlo Simulations 

The Monte Carlo simulations we carry out in this work have 
the advantage of computing both the continuum emission 



and the reflected iron line flux as measured by a distant 
observer. In this manner the equivalent width of the line can 
be estimated in a consistent manner. A simulation consists of 
sending N photons from the source isotropically with respect 
to its proper frame. iV c of these photons will reach directly 
the observer's collecting area and therefore contribute to 
the continuum flux. Ni of them will be received after having 
hit the disc and being reprocessed into iron line photons. 
These contribute to the observed flux of the line. For each 
photon received on the collecting area, the total spectrum 
is computed by considering the total power in a succession 
of freq uen cy bin s of extent Av. The following two sections 
( [2.4.1 and 2.4.2) describe how the total power in each bin is 
estimated, first for the continuum emission then the reflected 
emission. 



2.4-1 Continuum Emission 

Let us consider a photon emitted by the source with 4- 
momentum p s (equation ^ which reaches directly the col- 
lecting area at radius r with 4- momentum p . The redshift 
measured by the distant observer is therefore 



1 + z so = 



Ps ■ Vs 
Po ■ Vo' 



(23) 



Since v — 70 the denominator is simply the value of $ at 
the observer, "I> . The numerator is calculated using (JHj) and 
(§). We find 



1 + Z so 



A! 



L sin 9 sin P sin(</> p — . 



(24) 



where the subscript s indicates that the expression is evalu- 
ated at the source. 

We define W£ to be the continuum total power received 
in the collecting area of the distant observer in the range of 
frequency [v, v + Av\. We have 



(25) 



where the sum is for each received continuum photon, W Va 
is the received power per unit frequency attached to the 
photon, and the subscript o denotes quantities as measured 
by the distant observer. Expressing quantities on the right 
hand side of equation (^5|) in terms of quantities defined at 
the source, we obtain 



E 



\l + ZsoJ J v( 



0+Ai/)(1+ Z( , ) 



r+i j 
u. Ws 



(26) 



l + z so ) 



2.4-2 Reflected Iron Line 

We now consider a photon that leaves the source with a 4- 
momentum p s and hits the disc with a 4-momentum pd- The 
redshift measured by the co-rotating observer is 



1 + Z sd = 



Pd ■ v d 



(27) 



The numerator is given i n (B4h and the denominator can be 
calculated using (jpj) and (1131)7 We find 
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1 + z sd = 



A 3 — L sin 9 sin 9 P sin(0 p 



4\i 



cosh J7 — sinh U sin # p sin(0 p 



(28) 



where the subscripts s and d indicate that the expression 
is evaluated at the source and at the disc, respectively. A 
similar approach can be followed for the second part of the 
trajectory, from the disc to the collecting area of the distant 
observer. We find 



4\. 



1 + z do = 



cosh U 



(29) 



The radiation processes occurring on the disc are de- 
scribed in detail in George & Fabian (1991). Here, we are 
simplifying the situation significantly by ignoring the re- 
flected continuum emission so that only the reflected iron 
line is taken into account (see Section H for a discussion of 
the validity of this approach). We define Nl™(v d ,9i) to be 
the total number of photons impinging on the disc with en- 
ergy v d and incidence angle Qt per unit frequency and per 
unit time, as measured by a co-rotating observer (see Fig. pf). 
Following George & Fabian (1991), the number of fluores- 
cent photons per unit frequency and per unit time which are 
able to escape the disc is 



u d - "v d (yd,0i)g{0i)f(v d ), 



where, 

g(9i) = 



(6.5 - 5.6 cos 9i + 2.2 cos 2 9i) x 10" 
Vd - 1.: 



f(u d ) = 7.4 x 1(T 2 +2.5exp 



5.7 



(30) 

(31) 
(32) 



This analytical approximation is valid for v t < v d < u m , 
where u t = 7.1 keV is the energy threshold for trigerring 
the 6.4 keV iron fluorescent emission, and u m = 30 keV. 
As given in Ghisellini, Haardt & Matt (1994), the intensity 
of the emerging fluorescent emission is not isotropic and is 
here assumed to be proportional to cos# e ln(l + 1/ cos# e ), 
where 9 e is the outgoing inclination angle as measured in 
the co-rotating frame. Therefore iV°" follows the angular 
distribution 

N°f{e e )dO e = 2 cos (9 e In (l + —^—) sm6» e d6» e . (33) 
V cosfc» e / 

We now define Wl to be the iron line total power re- 
ceived in the collecting area of the distant observer in the 
range of freq uency [v, v + Au]. Following a similar reasoning 
as in Section 



2.4. lj we obtain 

2 



1 



1 + Z do 



N° 



(34) 



where v a = 6.4 keV, S u = 1 if v < v a /{l + Zdo) < v + Av, 
5 V — otherwise, and 



N° ut = g(0i) / f(v d )NT d dv d . 

J u t 

When expressed in the source proper frame, we have 

2 

X 

V 1 + Zrln ) 1 + Z*A " ' ' 

N, 



(35) 



^ V 1 + z do ) 1 + z sd 

\ 1 + z sd J 



(36) 



Source Height (r ) 

Figure 4. Ratio of the solid angle formed by the rays hitting 
the disc to the solid angle formed by rays escaping to infinity. 
The solid line is for a maximally rotating Kerr black hole while 
the dashed curve is for a Schwarzschild hole. The dotted-dashcd 
curve is for the ratio of solid angle formed by rays deflected to- 
wards the equatorial plane, including those lensed towards the 
event horizon. The dotted line indicates the ratio that would be 
obtained without the gravitational bending effect. 



3 SOURCE ON AXIS 

When discussing the fluorescent iron line profile and equiv- 
alent width, many authors have considered the situation 
where the primary X-ray source is located above the disc, on 
the axis of rotation (e.g. Martocchia et al. 1999, Matt et al. 
1992, Martocchia & Matt 1996; Reynolds & Begelman 1997, 
Young, Ross & Fabian 1998; Reynolds et al. 1999). This as- 
sumption has the advantage of simplifying both calculations 
and analysis since the axial symmetry of the problem is con- 
served. In this section we also follow this assumption, but it 
will be relaxed in Section ^j. 

3.1 Gravitational Lensing Effect 

As pointed out by Martocchia & Matt (1996), if the pri- 
mary source is very close to the hole, a significant fraction 
of the primary photon rays will be lensed towards the disc, 
decreasing considerably the direct flux relatively to the re- 
flected flux, as observed by a distant observer. In order to 
quantify this phenomena, the ratio of the solid angle formed 
by the rays hitting the disc to the solid angle formed by the 
rays escaping to infinity (f2disc/f2 osca p ) is plotted versus the 
height h — z s of the source in Fig. 0. The effect is negligi- 
ble in the case of a Schwarzschild black hole since a large 
amount of the lensed radiation is lost within the innermost 
stable orbit of the disc (i.e. 6 r g ). In the case of a maximally 
rotating Kerr black hole, where the inner radius of the disc 
can be as small as 1 r g , the effect becomes significant for a 
sufficiently close source, typically h < 6 r g as seen on Fig. ^. 
However, in the case of a source at height h — 2 r g , Q.disc is 
only a factor of ~ 5.5 bigger than fl csca p C . In this case we 
find that almost half of the rays bent towards the equatorial 
plane do not actually hit the disc but are strongly lensed to- 
wards the event horizon, resulting in a significant damping 
of the effect. 
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Figure 5. Predicted spectral line shapes for a central 
Schwarzschild black hole (L = 0), as seen by a distant observer at 
inclination angle 9q = 30 degrees. The primary source is located 
on the axis of rotation at height h from the origin. The solid, 
dashed and dotted lines are for h = 2.5, 5 and 30r 3 respectively. 
The area under each curve is proportional to the line equivalent 
width. 
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Figure 7. Predicted spectral line shapes for a maximally rotating 
Kerr black hole (L = 0.998A/), as seen by a distant observer at 
inclination angle 9q = 30 degrees. The primary source is located 
on the axis of rotation at height h from the origin. The solid, 
dashed and dotted lines are for h = 1.5, 5 and 30r 9 respectively. 
The area under each curve is proportional to the line equivalent 
width. 
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Figure 6. Same as Fig. H but for 6q = 60 degrees. 
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Figure 8. Same as Fig. [?] but for 0q = 60 degrees. 



3.2 Line Profiles 

Examples of line profiles obtained with our Monte-Carlo 
simulations for various values of h are given in Figs. ^| ^| 
and H for the maximum Kerr and the Schwarzschild cases. 
The predicted line shapes are for distant observers at incli- 
nation angles 9o — 30 degrees and 8o = 60 degrees. As 
already noticed in previous studies (e.g. Fabian et al. 1989, 
Dabrowski et al. 1997) , the line extends further towards the 
blue shifted part of the spectrum for larger inclination an- 
gles. Indeed, for large 9o the blue part of the line is boosted 
by the Doppler effect due to the rotation of the disc. On the 
other hand, the extension towards the red shifted part of the 
spectrum is controlled by the strength of the gravitational 
redshift which increases as the fluorescent emission happens 
closer to the central black hole. This is very clearly illus- 
trated here for the maximally rotating Kerr case (Figs [?], |^), 
where v/v a could be as small as 0.16 for a source located 
at h = 1.5 r g . Indeed, as the source gets closer to the hole 
it illuminates more intensively the inner part of the disc, 
resulting in strong gravitational redshift of the fluorescent 



line. In the Schwarzschild case as well (Figs, g, |6j) the red- 
dening of the line increases as the source gets closer to the 
hole. However, because the inner radius of the disc is 6 r g 
as opposed to 1.24 r g in the maximally rotating Kerr case, 
less strong reddening is predicted in the Schwarzschild case. 
For example, one can notice very little difference between 



the predicted lines for h 



2.5 r g . In partic- 



ular the red wing is not more extended in the h = 2.5 r g 
case because most of the illumination is focused within the 
minimum stable orbit. 



3.3 Equivalent Width 

The equivalent width (EW) of the iron fluorescent line is 
estimated as follows: 

rl. 



EW = 



Jo 



WB=u 



(37) 



Results are given in Figs. | and as a function of the height 
of the primary source in the maximally rotating Kerr case 
(L — 0.998M), in the Schwarzschild case (L — 0) and in an 
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Figure 9. Equivalent width of the iron fluorescent line as a func- 
tion of the height h of the source. The solid, dashed, and dotted 
lines are for L = 0.998M, L = 0.5M and L = respectively. The 
distant observer is at inclination 6>n = 30 degrees. 
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Figure 10. Same as Figjji] but for an observer at inclination 
8o = 60 degrees. 
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Figure 11. Predicted spectral line shapes for a maximally ro- 
tating black hole (L = 0.998A/), as seen by a distant observer at 
inclination angle O = 30 degrees and azimuthal angle <j> = n/2. 
The primary sources are located at p s = 2.5 r g and z s = 1 r g , in 
each case. The solid, dashed, dashed-dot and dotted lines are for 
<fts = 0, 4> B = n/2, 4> s = n and <f> s = 3n/2 respectively. The area 
under each curve is proportional to the line equivalent width. 

Fig.jl^, the gravitational lensing effect becomes more effec- 
tive when observed from a large inclination angle. Indeed, 
the line is strengthened by the lensing effect for both low 
and high inclination angle points of view. However, due to 
the same effect, the continuum emission received at large 
9 is strongly reduced compared to that received at smaller 
9 . For an observer at inclination 9 = 60 degrees, EWs as 
high as ~1.2 keV are predicted, which represents an EW en- 
hancement factor of about 8. This strong EW enhancement 
effect can also be seen in Fig. |^, where the area under each 
line is proportional to the line EW. However in the case of 
MCG-6-30-15 the disc is thought to have a inclination angle 
of about 30 degrees (e.g. Tanaka et al. 1995; Dabrowski et 
al. 1997). 



intermediate case where L = 0.5M. Fig. ^ is for a distant 
observer with inclination 9 = 30 degrees and Fig. |Io| is for 
#o = 60 degrees. We expect the strength of the line to be 
larger as the source gets closer to the hole since the pri- 
mary emission would preferentially be focused towards the 
plane of the disc. This effect is clearly visible in the case of 
a maximally rotating Kerr black hole (solid lines in FigsJ^ 
& |Io| ) . However, for lower values of the black hole spin and 
consequently higher values of the inner radius of the disc, 
this EW enhancement effect is lessened, as predicted in Sec- 
tion 3.1. In cases where the source is so close to the hole 



that most of the lensed illumination is lost within the inner 
radius of the disc, the effect is actually reversed and the EW 
diminishes as the primary source gets closer to the hole. In 
particular, no EW enhancement is predicted in the case of 
a Schwarzschild black hole. For a distant observer at incli- 
nation 9 — 30 degrees the maximum EW is predicted to be 
~290 eV, while the EW obtained for the classical case where 
no lensing effect comes into play (i.e. large h) is ~165 eV. 
Therefore, even for a primary source located just above the 
event horizon of a maximally rotating Kerr black hole, the 
EW enhancement is less than a factor of two. As seen in 



4 SOURCE OFF AXIS 

We have seen in Section ^ that the EW enhancement is 
limited by the fact that a significant proportion of the pri- 
mary illumination is 'lost' within the inner radius of the 
disc instead of contributing to the strengthen of the iron 
line. This limiting effect is particularly important when the 
primary source lies above the disc on the axis of rotation. 
Stronger EW enhancement should be obtained by relaxing 
this assumption and allowing the primary source not to lie 
right above the event horizon of the hole. In this section, we 
investigate this possibility in the maximally rotating Kerr 
solution (L = 0.998M). 

Since the axial symmetry of the system is here broken, 
we need to define the distant observer by both the inclination 
9 and the azimuthal angle <f>o- For simplicity, here we fix 
the observer at (j> — n/2 while allowing the source to be 
located at various azimuths <f> s above the disc, as suggested 
in Fig. |§|. Line spectra are given in Fig. [ll] as a function of (f> s 
for sources located at height z s — 1 r g above the plane of the 
disc and radius p s — 2.5 r g from the axis of rotation, close 
to the ergosphere radius at this latitude (1.93 r g ). Since the 
illumination takes place above the inner parts of the disc, 
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Figure 12. Equivalent width of the iron fluorescent line as a 
function of the height z s of the source. The primary source is 
at azimuth (f> s = 0, above the approaching side of the disc. The 
solid, dashed and dotted lines are for p s = 2.5, 3.5 and 5 r g 
respectively. The central black hole is a maximally rotating black 
hole (L = 0.998M) and the distant observer is at inclination 
O =3O degrees and azimuth <f> = tt. 
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Figure 13. Same as Fig. [12] but for an observer at inclination 
O =6O degrees. 

(p s = 2.5 r g , 4> a — 0, z s — 0.5 r g ), the observed EW can be 
as high as ~5.5 keV, which represents an EW enhancement 
factor of about 35. 



all these lines display strongly redshifted features. However, 
as one can notice on Fig. hll the shape of the spectral line 
varies significantly as a function of 4> a . For cf> s = the line 
is mostly redshifted and its flux peaks at around 4.5 keV, 
while for <f> s — -k the line is double peaked around 1.8 keV 
and 6.5 keV. In the former case the primary source is located 
above the approaching part of the disc where the emission 
is significantly boosted by the Doppler effect. In the latter 
case, the strongly redshifted red wing is the result of Doppler 
shift from the receding part of the disc, while the blue wing 
is due to rays that have reached the approaching side of the 
disc while spiralling down on to the disc. Lines for cf> s — tt/2 
and 4>s = Sn/2 are hybrids of the above two cases. 

Since the illumination occurs right above the Doppler 
boosted region of the disc for (j> s — 0, the predicted line EW 
is much larger in this geometrical situation than for other 
values of <f> a . This is illustrated in Fig. [y] where the spectral 
line integrals are proportional to their EW. For the rest of 
this section we assume that the primary source is always 
located at this privileged azimuthal angle <f> s — 0. Figs. [12 
and [H^ give the predicted iron line EW as a function of the 
primary source height z s and radius p s . In the case of an 
observer at inclination 8 o =30 degrees, the observed EW can 
be much larger than for a source located on the axis of ro- 
tation. For example, as seen on Figjl^, for a primary source 
located at (p a = 2.5 r g , <f> s — 0, z a — 1 r g ) the predicted EW 
is ~1 keV and if the source is as close as z 3 = 0.5 r g to the 
disc the EW can be as high as ~1.5 keV. In this latter case, 
the EW is enhanced by a factor of ~10 of its classical value. 
However, in order to reach such high EW the primary source 
needs to be located above the very inner parts of the disc. 
For example, in the case of a source located at (p s — 5 r g , 
4> s — 0, z s — 0.5 r g ) the predicted line EW is only ~500 eV. 
As a consequence, very little EW enhancement is expected 
in the case of slowly rotating black holes. As expected, for an 
observer at inclination 6> o =60 degrees the predicted line EW 
is much larger than for 6* o =30 degrees. As seen on Fig,|l3|, 
in the extreme case where the primary source is located at 



5 DISCUSSION 

In Section ^| we have presented a new fully relativistic ap- 
proach to calculating photon paths and energies in the 
Boyer-Lindquist metric. The geodesic equations obtained 
in this approach are simple first order ordinary differen- 
tial equations which allow easy and stable numerical in- 
tegration. Furthermore, the way the photon 4-momentum 
is parameterised (equation allows us to be very clear 
regarding observable quantities such as the frequency and 
incoming/outgoing angles of the light ray, as measured by 
an observer. We believe that our theoretical approach is an 
improvement compared to more standard methods usually 
employed in similar applications (e.g. Reynolds et al. 1999). 

We found in Sections ^ and ji] that large values of the 
observed iron line equivalent width (EW) may be explained 
by general relativistic effects when the primary X-ray source 
lies above the very inner part of an accretion disc (within 
~5 gravitational radii) around a maximally rotating central 
black hole. When the primary source is located on the axis 
of rotation, EWs up to ~300 eV only can be accounted for 
by this model for a disc at inclination 30 degrees. This is a 
much lower prediction than that of the early work of Mar- 
tocchia & Matt (1996) who found EWs up to ~1500 eV for 
a comparable case. Our results are however in good agree- 
ment with the very recent work of Martocchia et al. (1999). 
At inclination 60 degrees EWs up to ~1.2 keV are predicted. 

By allowing the primary source to be located off the axis 
of rotation, we found that much larger values of the EW are 
obtained when the source is located above the approaching 
side of the disc (up to ~1.5 keV and ~5.5 keV for inclina- 
tion 30 and 60 degrees respectively). However, in the most 
extreme cases, a reflected continuum component should also 
be included in the model, in addition to the reflected iron 
line (Lightman & White 1988; Pounds et al. 1990; Matt et al. 
1991; George & Fabian 1991 ). When the reflected emission 
is not enhanced by the lensing of the primary light rays, 
the reflected continuum accounts for about 10 percent of 
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the observed continuum flux at 6.4 keV. One would there- 
fore expect the reflected continuum to dominate when the 
enhancement factor reaches ~10. Beyond this value, the con- 
tinuum flux observed at 6.4 keV increases together with the 
iron line flux and therefore the line should no longer gain in 
strength. As a result one should regard values of Fig. [l3] as 
indicative of the strength of the lensing enhancement effect, 
while the actual observed iron line EW should not be larger 
than ~ 1600 eV (i.e. 10 times its classical value). 

Results found in Section ^ may be used as arguments to- 
wards explaining the large iron line EW observed in Seyfert- 
1 galaxies. Such an explanation may be valid for the EW 
enhancement observed in the case of MCG-6-30-15 where it 
is likely that strong emission processes take place very close 
to the hole, at least during specific periods (see Iwasawa 
et al. 1996, 1999). However, as highlighted by Reynolds & 
Fabian (1997), in the case of more typical objects it is less 
probable that strong emission would occur so close to the 
hole. Reynolds & Fabian (1997) found that the relative mo- 
tion between the disc and the source may significantly affect 
the EW of the line, simply by special relativistic arguments. 
This could allow EW enhancement even when the emission 
occurs further away from the hole. Such an effect could be 
investigated in our model by allowing the source to have a 
peculiar velocity. For example, a simple model where the 
source moves with a constant angular velocity could be as- 
sumed (Yu & Lu 1999). 
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APPENDIX A: GENERAL RELATIVISTIC 
APPROACH 

This work has been carried out using the gauge-theoretic 
formalism described in detail in Lasenby et al. (1998). How- 
ever, the theoretical calculations presented in Section ^ and 
Appendix |§| could in principle be carried out following a 
more standard general relativistic approach by making use 
of orthogonal tetrads (see e.g. Misner et al. 1973). A suitable 
tetrad frame in the Boyer-Lindquist metric (equation ^) can 
be defined as follows: 



Dotani T. 
Makino F.. 



= 

W = 



i ) * 

0, ,0,0 

p 



L 



0,0, -0 

p 

Lr sin 2 9 
P 



0,0,- 



(Al) 
(A2) 
(A3) 
(A4) 



where p = ... 3 are the usual coordinate indexes, corre- 
sponding here to the Boyer-Lindquist coordinates (t, r, 9, 
</>). This tetrad can be identified with the coordinate frame 
{e M } defined by equations ^ to (^). The orthonormal tetrad 
corresponding to the {7 M } frame (equation ^) is given by 



(7o)" 



■ 2 + L 2 

pA5 



,0,0, 



L 



(A5) 
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-L sin 8 sin d> A'sin6>cosrA 



(71)" 



(72)" 



P 

cos 6 cos ( 



p ' psin# 
L sin 6 cos 6 A 2 sin sine 



cos # sin c 



p psm( 



(A6) 



(A7) 



(73)" = (o ) ^^,^£,oj. (A8) 

The equations given in Section ^ can now be translated into 
this general relativistic context. For example equation Q of 
Section 



2.1 



should be read as follows 

P M = $(70)" +$ sin e p cos <?!>p (71 )" + 

$ sin P sin P (72)" + $ cos 6> p (73)" . (A9) 

While the photon geodesic equations t, r, 8, (j> given in Ap- 
pendix |b| are obtained by the following equation 

p dX 
where 



(A10) 



(dx)>* = (dt, dr, dO, &<t>) . (All) 
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APPENDIX B: PHOTON GEODESIC EQUATIONS 



When solving the geodesic equations in the Boyer-Lindquist metric for the photon 4-momentum given in Section 2.1 (see 
Lasenby et al. 1998), we find 



dA 

dr 
dA 
d0_ 
dA 



$A? 



L 2 + r 2 + LA 2 sin 9 sin 9 P sin(</> p 



cos 9 P cos 9 + sin 9 sin 9 P cos(0 p 



(Bl) 



cos V sin Up cos 



sm cos ( 



i 



dA pA*sin6» 
and for the photon 4-momentum 



A' sin 9 P sm(4>p — </>)+ L sin 9 




d£ 
dA 



dA 

dffp 
dA 

where 



A = 2L sin # p sin(<j!> p 



cos 9 cos 9 P + sin 9 sin # p cos(<^ p — 0) 



+ LA 



2 2 

r 



sin 



(2p ; 
b) sin 20 



2 tan 9 V 



r A 2 + B) ~ L 2 cos 6> cos ( 



Lp'A" 



(B2) 



S[tan( 



) cos # p — sin # p ] + LA sin 9 sin(0 p — 0) 



B = Ar + p 2 A~ 



4 ( 



)sin6> + 2A 2 , 
M - r - A' 



We note the very simple form of the ordinary first order differential equations obtained for t, r, 9 and 4>, where overdots 
denote differentiation with respect to A. The $, 9 P and (b p expressions are slightly more complex but still are first order ordi- 
nary differential equations. This approach is therefore useful because of the computational advantage of first order equations 
in terms of efficiency and stability. 



An interesting aspect of the geodesic equations (Bl) and 



is that they do not rely on conserved quantities such 



as the energy, angular momentum or Carter constant along the photon path. As a consequence, one can actually compute 
these quantities and verify that they are indeed conserved along the trajectory. This allows one to check easily the analytical 
expressions (Bl) and (B2) as well as the corresponding numerical code. The conserved energy E p and angular momentum J p 
of a photon geodesic are evaluated by making use of the Killing vectors K t and K$ associated with the temporal and axial 
symmetries of ([[]) respectively. We have 



K t = 
and 



A 2 



L 

— f 

pr 



LA? sin 2 9 



-e t + 



r 2 + L 2 



pr 



-e<j,. 



Therefore 
E p =p-Kt = $ 
and 

J p — p ■ K$ — — <£> sin i 



L sin 9 sin 9 P sin(0 p — (f>) + A' 



(L 2 + 



r smS„sm 



+ LA 2 sin( 



(B3) 



(B4) 



(B5) 



(B6) 



The derivatives E p and J p with respect to the affine parameter A can be obtained. As expected, substituting for the expressions 
of r, 9, cj>, $, 9 P and (b p found above leads to 

This paper has been produced using the Royal Astronomical Society/Blackwell Science M]gX style file. 
© 1997 RAS, MNRAS 000, 



